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The joint analysis of galaxy-galaxy leasing and galaxy clustering is a promising method for in¬ 
ferring the growth function of large scale structure. This analysis will be carried out on data from 
the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the 
tangential shears of background galaxies induced by these foreground lenses. We develop a practical 
approach to modeling the assumptions and systematic effects affecting small scale lensing, which 
provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize 
the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement 
errors, we study how external priors on different subsets of these parameters affect our growth con¬ 
straints. Degeneracies within the HOD model, as well as between the HOD and the growth function, 
are identified as the dominant source of complication, with other systematic effects sub-dominant. 

The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of 
the galaxy sample that we employ. We conclude that DES data will provide powerful constraints 
on the evolution of structure growth in the universe, conservatively/optimistically constraining the 
growth function to 7.9%/4.8% with its first-year data that covered over 1000 square degrees, and to 
3.9%/2.3% with its full five-year data that will survey 5000 square degrees, including both statistical 
and systematic uncertainties. 
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f Corresponding author lise@slac.stanford.edu 



3 


I. INTRODUCTION 


Evidence from multiple probes now points to an accelerated expansion of the Universe. Distant Type la supernovae 
are fainter than they would be if the Universe were decelerating mm-: patterns in the anisotropy of the Cosmic 
Microwave Background (CMB) have long been consistent with acceleration and now offer solid independent evidence 
[3]; the scale of Baryonic Acoustic Oscillations in the late-time galaxy distributions also points to acceleration [3]. 
Other measurements, while not providing stand-alone evidence, are nonetheless consistent with the notion that the 
deceleration predicted by Einstein’s theory of General Relativity without a cosmological constant is not occurring 
today. For example, measurements of growth of structure using the abundance of massive clusters of galaxies [Si, 
as well as weak lensing [7] have been found to be consistent with a model in which dark energy driving acceleration 
contributes roughly 70% of the energy density of the Universe. The physical nature of the mechanism driving this 
accelerated expansion, however, is still to be determined. 

A major goal of the Dark Energy Survey (DES) is to understand that mechanism by measuring the growth of 
large scale structure. Because different models predict distinct histories of structure growth in the late-time universe, 
constraints on growth history can lead to constraints on the mechanism responsible for cosmic acceleration. We 
expect the most precise constraints to be obtained using combinations of several probes (e.g., see 0), which increase 
the overall signal-to-noise and break parameter degeneracies - both among the cosmological parameters of interest 
and the nuisance parameters that quantify systematic effects. The example we focus on here is the combination of 
measurements of galaxy-galaxy lensing and clustering of the lens galaxy sample, which has been suggested in the past 
few years by liiin]. By constraining the growth function with such a combined analysis, we not only constrain the 
parameters of the “standard model” of cosmology but also can detect possible deviations from the robust predictions 
of General Relativity and smooth dark energy models. 

In particular we implement the approach proposed in Yoo and Seljak [9], which combines small-scale galaxy-galaxy 
lensing with large-scale clustering, on mock data sets designed to resemble that from DES. On large spatial scales, the 
galaxy overdensity is proportional to the overdensity in the total matter distribution, with the relation between the 
two over-densities captured by a single number, the linear bias parameter (e.g., see m) which is related to the masses 
of halos hosting the galaxy sample. On small spatial scales the relation between galaxy and dark matter distribution 
is non-linear. The small-scale dark matter distribution is assumed to follow that of a spherical halo with a universal 
mass profile, and the distribution of galaxies within a halo is commonly described by Halo Occupation Distributions 
(HOD) [12]. HODs are used extensively to model galaxy-galaxy lensing and clustering [e.g. [101 [E], and have been 
successfully applied in recent joint analyses of galaxy-galaxy lensing and clustering [IMISj- The insight of |5] was 
that one could apply a step-by-step method to address these different scales and corresponding physics, starting by 
fitting lensing on small scales with a simple mass profile. Then, the inferred mass could be used to understand the 
large scale bias of the lensing galaxies, turning large-scale galaxy clustering measurements into direct probes of the 
underlying clustering of matter. By carrying out this two-step analysis with lensing galaxies in multiple redshift bins, 
one might therefore be able to measure the growth of the structure. 

Here, we implement this idea anticipating a near future application to DES data. We find that the simple two-step 
approach needs to be tweaked, and that parameters accounting for the HOD modeling as well as for systematic 
effects need to be introduced into a full one-step analysis that includes both sets of measurements - clustering and 
lensing - in one data vector. We develop a full analysis pipeline for this method and forecast its constraining power 
at different stages of the Dark Energy Survey. We employ a joint model for key systematic effects such as halo 
model assumptions and photometric redshift errors, allowing for both probes to contribute information on, and best 
constrain, the underlying assumptions and parameters under a realistic setting. In addition, to correctly account for 
the correlation between probes, we utilize the full joint covariance matrix of the two probes. We test and validate 
this pipeline in simulated data designed to closely mimic that obtained and expected from DES. This pipeline will 
be applied to DES data, so one of our goals here is to test the pipeline and the underlying algorithm: how can we 
optimize this joint analysis on actual survey data given the statistical uncertainties and likely sources of systematic 
error? Which systematic effects are most important to model accurately and which do not affect the final cosmological 
constraints? Most generally, how accurately should we expect to be able to extract information about the growth of 
cosmic structure? 

The plan of the paper is as follows. Section H contains a description of the implementation: the halo model based 
formalism and the choice of our parameter set. The mock catalogs, measurements and tests are presented in Section 
HI. In Section IV we describe our likelihood analysis and details on model parametrization. We present and discuss 
our results in Sections V and VI. 
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II. MODELING 
A. Motivation 

The focus of this paper is to develop a pipeline that will extract information about the growth function from 
small-scale DES galaxy-galaxy lensing and large-scale DES galaxy clustering. 

A first attempt to implement the method described in Ref. [S] would be to: 

1. Select a galaxy sample with a given luminosity cut, with a parametrized model for the mass-luminosity relation 
and redshift range. 

2. Eit the halo-scale galaxy-galaxy lensing data, 7t(0), with a halo mass profile to extract an estimate of the mean 
mass of the sample. 

3. Determine the large-scale halo bias for that mass using fits from numerical simulations, e.g. m- 

4. Measure the angular correlation function, w(0), of the galaxy sample. 

5. Using the inferred halo bias and external priors from, e.g., Planck [18) . simultaneously ht the correlation function 
to a set of cosmological parameters including the growth function. 

Writing these steps down immediately reveals a number of problems. In order to carry out each of Steps 1-3, a 
distance-redshift relation is needed, which depends on cosmology. So in principle, one cannot fix this relation and 
then at the final step fit for cosmological parameters. Second, redshift bins will be determined using DES colors so 
will be subject to photometric redshift errors, and these affect the fits in Steps 1, 3 and 5. Therefore uncertainties in 
photometric redshifts must be treated simultaneously. Finally, some information is needed about the mass-luminosity 
relation and particularly about the fraction of galaxies that are satellites instead of central galaxies. For these purposes 
a more sophisticated analysis is needed even at the outset. 

We aim to maintain the basic idea of Yoo and Seljak [9] of combining small-scale galaxy galaxy lensing with 
large-scale galaxy clustering, while addressing the above issues. Our starting point then is the joint data vector that 
includes both 7 t( 0 ) and w(0) for the luminosity-threshold galaxy sample. To extract predictions for these statistics, we 
employ a halo model |19j in combination with HOD modeling. Specifically, we dehne halos as spherical overdensities 
of Am = p!P m = 200, and assume their densities follow the Navarro, Frenk, & White (NEW) profiles [20] with the 
Duffy et al. jH] mass-concentration relation. We use Tinker et al. iniiii] fitting functions for the halo mass function 
and halo mass-bias relation, respectively. We then jointly model both w{9) and 7 t( 6 *) from this halo model picture, 
with added ingredients for systematic effects such as photometric redshift errors and multiplicative shear calibration. 
In addition to the parameters associated with the HOD modeling, the set of systematic effects, and cosmology, we 
introduce growth scaling parameters, denoted A^, to freely scale the amplitude of the growth function in each redshift 
bin, rendering our analysis capable of both constraining the growth function and detecting potential deviations from 
ACDM structure growth. A key ingredient of this analysis is the full joint covariance matrix of the joint data vector. 
In treating the joint likelihood, the full joint covariance matrix allows for a proper accounting of the information in 
the joint data vector, especially with its off-diagonal blocks representing covariances between the two probes. 


B. Halo Occupation Distribution 

When we measure the tangential shear induced by stacked foreground halos, what is the best way to characterize 
our sample? Simply fitting for a single value, i.e. the mean halo mass, is not optimal because it does not fully 
represent the underlying mass distribution of halos, thereby leaving out information. Rather, directly modeling that 
underlying mass distribution by means of a halo mass function will yield a more realistic characterization of the 
sample. Furthermore, we observe galaxies, not halos, so in addition to the mass function we also need a recipe that 
connects galaxies to halos. Going from a halo mass function to a galaxy distribution requires an HOD model that 
describes the relation between galaxies and halo mass in terms of the probability P{N\M\^) that a halo of given mass 
Mh contains N galaxies. We separate galaxies into central and satellite galaxies. By definition, a halo contains either 
zero or one central galaxy, and it can only host satellite galaxies if it contains a central galaxy, which motivates the 
form |23| 


(?V(Mh)) = (iVc(Mh)) (1 + (iVs(Mh))), 


with the average number of central/satellite galaxies in a halo of mass M-^. 


(1) 
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FIG. 1. An example of the average number of central/satellite galaxies, (Ac/s(Afh)), calculated from Eq. [^with parameter 
settings logMmin = 12.36, log Mi = 13.69, friogM = 0.32, a = 1.28. These parameter values are selected to match our hducial 
default values presented in Table [I] The dashed and dotted black lines respectively represent the central and satellite galaxy 
counts, with the solid black line showing their sum, i.e. the total number of galaxies in a halo of mass Mh. The solid and 
dashed red lines respectively represent the satellite and total counts using log Mo = 8.35 in addition, i.e. Eq. [^before our 
simplihcation. For the galax;y sample under consideration, the effect of the satellite cut-off mass scale Mq is negligible. 


For a luminosity-threshold sample (with absolute r-band magnitude Air 
is commonly parameterized as [e.g.. 123) 


< Ail), the HOD for centrals and satellites 


(7V,(Mh|M‘)) 


1 [l , r log Mh - log Mmin 

^ \ ^\ogM 

/Mh-Mo\“ 

{ M[ ) ’ 


( 2 ) 


with model parameters Almin, M(, uiogM, ot, and all mass parameters in units of Mq/H. The central galaxy occupation 
function is a softened step function with transition mass scale Almin, which is the halo mass in which the median 
central galaxy luminosity corresponds to the luminosity threshold, and softening parameter uiogM which is related to 
the scatter between galaxy luminosity and halo mass. The normalization of the satellite occupation function, M{, and 
cut-off scale Mq are related to Mi, the mass scale at which a halo hosts at least one satellite galaxy {{Ns{Mi)) = 1)); 
finally a is the high-mass-end slope of the satellite occupation function. This parametrization was found to reproduce 
the clustering of SDSS [24] and CFHTLS [25] galaxies well over a large range of luminosity thresholds and redshifts. 
To simplify this model and reduce the number of fit parameters, we ignore the satellite cut-off scale Mg = 0 and 
use a four-parameter model for luminosity threshold samples. Fig. [^ illustrates our HOD model, exhibiting the 
soft low-mass threshold determined by logM^in and tJiogM, the satellite onset dictated by log Mi, and the rapid 
increase of satellite counts at the high-mass end governed by a. Section P^I Ej describes the halo model that relates the 
HOD to galaxy-galaxy leasing and clustering observables. Throughout we use the mean value of central and satellite 
occupation, and assume that satellite galaxies are Poisson distributed both in number and in position with respect to 
halo centers to calculate second moments. 


C. Photometric Redshift Uncertainties and Shear Calibration 

The redshift distribution of galaxies plays a key role in projecting the 3D information to the 2D observables w{Q) 
and 7t(0), as well as in interpreting the tangential shear of a source galaxy image by a lens galaxy. In photometric 
surveys like DES the true redshifts of observed galaxies are not available; instead, redshift values are estimated from 
a galaxy’s brightness in different colors, known as photometric redshifts, or photo-z’s, Zph. Galaxies with photometric 
redshifts within a given range are lumped into a photometric redshift bin. To infer the true redshift distribution of 
this bin, we convolve the conditional probability function p(z\zph) with the photometric redshift distribution n{zph) 
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to calculate the true redshift distribution of the i-th photometric redshift bin ni(z): 


ni{z) = 



dzphp(2rphk)n(2;ph) ■ 


(3) 


We assume a Gaussian distribution of photo-zs around a true redshift value with the redshift-dependent standard 
error ct = tTz(l + z) and constant offset [e.g.. 1^. 


P{zph\z) 


,— exp 


{z - Zph - 

2 ct 2 


(4) 


This model is a somewhat idealized picture, as in reality complex galaxy spectra give rise to complicated, non-Gaussian 
photo-z distributions. Here we choose this simple parametrization to manifest the most important modes of error in 
photo-z’s, also noting that the Gaussian assumption holds well for our expected candidates for lens galaxies, namely 
Luminous Red Galaxies (LRG’s). 

In addition, we consider a multiplicative calibration of the observed tangential shear as a potential source of 
systematic effects, so that the true shear is related to the observed shear via 


(7t(6'))true = (1 + (7t(^))obs ■ 


(5) 


D. Growth Function Scaling 


At the linear level, the growth of structure in the universe is described by the growth function D[z), normalized to 
be unity at z = 0. For example, in terms of D(z), the matter power spectrum P{k, z) is 


P{k,z) = D^{z)P{k,0), 


( 6 ) 


which then enters various structure-related quantities such as the variance of matter density fluctuations on a scale 
R, crfl(z), and subsequently the mass function dn/dM\y. For a standard flat LGDM cosmology, D{z) is given by 


^ACDM(^) 


H{z) 

/■~ dz'(l + z') 

■ /■°° dz"(l + z"y 

Ho 

A H^z') 

[Jo H^z") \ 


(7) 


1 /2 

where H{z) = Hq (f2M(l + z)^ + Ha) with the present-day Hubble constant Hq and density parameters Hm, Ha. 
Therefore, these three parameters uniquely define the growth function in the LGDM scenario. Therefore, in order to 
capture sensitivity to possible anomalies in the growth function, we introduce free scaling parameters Ai defined 

b,{z)=A,D^^^^{z), ( 8 ) 


that scales the growth function for the Tth redshift bin in our galaxy sample. The ensuing constraints on A^ capture 
the sensitivity of the combined probes to the amplitude of fluctuations at the redshift of interest. If the Ai are fonud 
to differ from one at a significant level, then LGDM would be ruled out. More generally, modified gravity models 
make different predictions for growth than do dark energy models, so independent measures of growth such as the Ai 
are extremely valuable ways to distinguish between these competing ideas for the cause of the cosmic acceleration. 


E. Observables 


1. Large-Scale Galaxy Clustering 


Our analysis uses the two-point function of the galaxy distribution on scales larger than individual halos. The 
angular power spectrum of galaxies in a given redshift bin i then depends on the linear matter power spectrum via 

Cl^il) = J dzH{z)x-\z)Wl,{z)P{k = l/x, z), (9) 

where x(z) is the comoving distance out to redshift z, and the galaxy window function Wgg{z) in bin i is 

WgAb = ( 10 ) 

Tii 
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with ni{z) the redshift distribution inferred from photometric estimates (see Eq. [^, and normalization factor rii = 
/ dzni{z). The mean galaxy bias bg(z) is given by 

1 ^OO J 

6g(z) = ^ / dM^ (iV(Mh|X)) (11) 

riM Jo dMi, 

with X = {Mmin, (TiogM, a} representing the HOD parameters defined in Eq. Here, dn/dM'o and b\y{M\J) are 
the halo mass function and the halo mass-bias relation from Tinker et al. [25] and Tinker et al. HZ], respectively. Note 
that as these quantities depend on cr(i?, z), they are affected by the growth scaling parameters Ai. The normalization 
parameter fiM is given by 


pOO J 

nM= / dMh^(iV(Mh|X)). ( 12 ) 

Jo 

In the flat sky limit, our observable w{9) is related to Cgg{l) as 

^(0) = j ^-^Cgg{l)J0{ie), (13) 

where Jo{W) is the zeroth-order Bessel function. 


2. Small-Scale Galaxy-Galaxy Sensing 

The measured tangential shear {"HJ {9)) of foreground galaxies in redshift bin i and source galaxies in redshift bin 
j is related to the Fourier transform of the tomographic galaxy-convergence angular power spectrum, Cg^^iJ) by 


{i?m = /^c'^(o^ 2 (z 0 ), 


(14) 


with J 2 the second-order Bessel function. 

The angular galaxy-convergence power spectrum is an integral over the 3D galaxy-mass power spectrum; in the 
small angle Limber approximation. 


Cli{l) = I dzx-\z)^^Wi{z)PgUk = l/xz). 

J 


Here, the lensing window function (z) for source bin j is 

TTi(z) = 


Pm{z) 


(l + ^)ScritW’ 


where the critical surface density (z) of source bin j is given by 


(s;.,.) 


AirGxiz) 


+ z 


i-x(z) 


XiZs)/\ ’ 


(15) 


(16) 


(17) 


with (x~^('2s)) the mean inverse comoving distance to the source galaxies in source bin j. 

It remains to compute the 3D galaxy-mass spectrum, which we describe using the halo model and HOD. For 
this analysis we will ignore the contribution of sub-halos and model the lensing signal around satellite galaxies with 
mis-centered NEW halos. Since we focus on small scales, we consider only the one-halo term: 


1 /* d 

P^!i,{k,X) = P,^{k,X) + P,^{k,X) = / dMi,Mi,Ui,{k,Mi,)^ [(iVc(Mh|X)) + {N,{Mi,\X)) u,{k, Mi,)] , 

PmUM J aMh 

(18) 

where Uh{k, Mi,) is the Fourier transform of the halo density profile of mass Mh, and Us{k, Mi,) the Fourier transform 
of the spatial distribution of satellite galaxies within the halo. Here, we assume that the distribution of satellite 
galaxies follows the NFW profile by letting Ug = Uh, and also that central galaxies are located at the exact halo 
centers, i.e. without mis-centering. 














III. MOCK DATA 


A. DES Data Stages 

DES is an ongoing wide field multi-color imaging survey that will cover nearly 5000 square degrees of the southern 
sky with a limiting i-band magnitude of 24 by Spring 2018. Its images come from the Dark Energy Camera 
a 3 square degree imager on the Blanco Telescope near La Serena, Chile. Images taken by the camera to roughly 
comparable depths will be obtained in g, r, i, z, Y bands, which will be used to characterize the positions, redshifts, 
and shapes of about 300 million galaxies. Pre-survey science verification data was taken from December 2012 to 
February 2013 and processed in Fall 2013. This data set, named the SVAl data release, covers about 150 square 
degrees to a limiting magnitude rrir 24 in the r band. The first year of science observations, referred to as Yl, 
has been released, covering over 1000 square degrees to roughly 0.5 magnitudes shallower depth. The complete DES 
dataset, to be achieved with 5 years of full data taking, is referred to as the Y5 dataset. 

To test our modeling discussed in Section|nj we construct a number of fiducial datasets from numerical simulations. 
In this work, we consider two different DES data stages, namely the DES Yl and DES Y5 stages. The full range of 
our pipeline is tested using simulated likelihood analyses with DES Y 1-like and Y5-like survey parameters and mock 
covariances from the DES Blind Cosmology Challenge (BCC) simulation results. 


B. Mock Survey Setup 

We make use of the DES BCC mock galaxy catalogs developed for the DES collaboration (Busha et al. [28]) to 
construct our mock surveys. From these DES galaxy mock catalogs, we construct luminosity-thresholded lens galaxy 
samples over two redshift bins, namely at 0.3 < z < 0.4 with Air < —21.5 and and at 0.4 < z < 0.5 with Air < —22.0. 
In order to obtain a realistic galaxy-galaxy leasing signal from mock catalogs with finite mass resolution the host halo 
of the lens galaxy needs to be resolved. Hence the luminosity thresholds for our lens samples are chosen such that 
central galaxies are located in resolved halos. 

The source sample is selected from the DES shear mock catalog by additionally imposing mi < 23.0 for the Yl 
source sample, and mi < 23.5 for the Y5 source sample to model the different depth of these two survey stages (Huan 
Lin, private communication). The resulting source catalog has an effective source density of 4.34 galaxies/arcmin^ 
(2.70 galaxies/arcmin^) for Y5 (Yl). While the mi < 23.5 is shallower than the nominal survey depth of Wi ^ 24.0, 
the resulting source galaxy density for Y5 is comparable to that of current shear catalogs for the SVA data |5^. We 
divide these background sources into three source redshift bins, 0.5 < z^ < 0.8, 0.8 < z^ < 1.1, and 1.1 < zf < 2.0. 
Figure shows the resulting redshift distributions of lens and source galaxies. The source tomography bins contain 
^gai ~ {1.25,0.46,0.28} galaxies/arcmin^ for our DES Y5 model, and = {0.65,0.18,0.11} galaxies/arcmin^ for 
our DES Yl model. 




FIG. 2. The lens redshift distributions (left) and source redshift distributions (right) are presented. For the lenses, we show 
the measured redshift distribution (dashed) and the deduced true redshift distribution (solid) for the two lens bins of our mock 
catalog. For the sources, we show the Yl and Y5 redshift distributions, with color-filled regions indicating the three source 
bins used. 
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C. Measurement Vector 

For our simulated likelihood analysis, we generate a measurement vector from our modeling framework assuming 
a set of fiducial default values for parameters. That is, we use the output of our prediction codes for 7 t( 6 >) and 
w{9) under the hducial default parameter settings as our measurement vector. We will refer to this as the simulated 
measurement vector. This, by construction, ensures that we can examine the information content of the proposed 
method, which is the goal of this paper, independent of discrepancies between simulations and theoretical models. 

We choose the small-scale lensing data vector to range from 1 to 6 arcminutes across 9 logarithmic bins, and the 
large-scale clustering data vector to range from 15 to 150 arcminutes across 10 logarithmic bins. This choice of scales 
separates the one-halo regime from the large-scale, linear clustering regime, and cuts out the transition and weakly 
non-linear clustering regimes, where the theoretical modeling uncertainty is the largest. 


D. Covariance Estimation 


We approximate the survey geometry of the Y1 and Y5 DES footprint as rectangles of 1000 and 5000 square degrees, 
respectively. We use the tree code treecor [SD] to calculate 7t(0) and w{9) (using the Landy-Szalay estimator [3T] 
with uniform random mocks for the latter), and measure the joint covariances by the bootstrap-with-oversampling 
method of Norberg et al. |32] . using 20 square degree patches and an oversampling factor of 3, yielding 

1 ^ - 

Cov(d*, dj) = (df - d,) {d) - di) (19) 

k^l 

with the joint data vector d = z^)), the k-th bootstrap 

realization, N = SlVpatch the number of bootstrap samples, and d the mean data vector calculated as 


d = 


N 


N 






( 20 ) 


We estimate the joint clustering and galaxy-galaxy lensing covariance for the two lens bins separately, and assume a 
block-diagonal total covariance matrix for the combination of multiple lens bins. 



FIG. 3. Left: Correlation matrix of galaxy clustering and galaxy-galaxy lensing with a single source bin for illustrative purposes. 
A large range of scales (1' — 100') is shown with a larger lens sample to reduce shot/shape noise and highlight the correlations 
of density modes. The black box indicates the range of scales considered in this analysis. Right: Correlation matrix of galaxy 
clustering and tom ograph ic galaxy-galaxy lensing for the DES Y5 0.3 < z < 0.4 lens sample and range of scales considered 
in this analysis (cf. Ill Cl. The panels marked as Zj correspond to the tangential shear measurement vector at the jth source 
redshift bin. 









































10 


The left panel in Figure shows the correlation matrix of clustering and galaxy-galaxy lensing over a larger range 
of scales than those considered in this analysis to illustrate the correlation between scales and probes. The black box 
indicates the range of scales considered in the analysis. Note that this covariance matrix is based on a larger lens 
sample in order to reduce statistical noise and highlight the underlying correlations due the correlation of density 
modes and the fact that galaxy-galaxy lensing and clustering both probe the underlying matter density field. In the 
right panel, we show the actual correlation matrix for our Y5 data vector in a single lens bin, with the tangential 
shear measurements from the three source bins marked as Zi. We observe reduced off-diagonal covariances, as shape 
noise and shot noise (respectively) are higher for the tomographic galaxy-galaxy lensing and the clustering of the lens 
galaxy sample used in our analysis. 

Note that while we choose the mock lens galaxy samples used in the covariance estimation to be similar in mass 
range and number density to the fiducial lens galaxy sample used for generating the measurement vector, the match 
is not exact. In order to adjust for the difference in signal strength to leading order, we rescale the clustering auto¬ 
covariance, clustering - galaxy-galaxy lensing cross-covariance, and galaxy-galaxy lensing auto-covariance by their 
respective scaling with galaxy bias, i.e. by ( 6 fid/&mock )^5 (^fid/^mock)^, and ( 6 fid/^mock)^ respectively. Here, 6 fid is the 
galaxy bias calculated for the synthetic measurement vector, and 6 inock is the galaxy bias measured from the mock 
data. This covariance rescaling is equivalent to performing an analysis using the original covariance with rescaled 
HOD-derived data vector (w, 7 t) —>■ ((^mock/&fid)^w^, (^mock/^'fid) 7 t), and does not change the shot noise level. The 
latter is difficult to adjust in real-space covariances as (due to a mixed cosmic variance and shot noise term) it affects 
all covariance elements differently [33]. Since the number densities corresponding to our fiducial HOD parameters for 
the lens sample are higher than the number densities of the mock samples, this is a conservative rescaling and may 
overestimate statistical errors. 


IV. LIKELIHOOD ANALYSIS 


A. Overview 


we perform a Markov Chain Monte Carlo 
(MCMC) likelihood analysis to forecast how well this analysis can constrain model parameters under various data 
stages of DES. To generate the simulated measurement vector used in these analyses, we must define a set of fiducial 
default values for model parameters. These fiducial defaults represent our best-guess estimates for the model parame¬ 
ters characterizing actual DES data. In addition, for the likelihood analysis, a set of priors for these parameters must 
be assumed. Priors allow us to include information outside of our pipeline, either from DES or from external results, 
that strengthens our constraining power. It is important to note that our priors on the HOD and systematic effects 
parameters represent the constraining power we expect to obtain with DES data outside of our pipeline, even when 
we benchmark our estimates from external results where we lack existing analyses of DES data. Also, since we use a 
simulated measurement vector without random errors here, i.e. without introducing further random fluctuations to 
the output from the prediction code, we focus on investigating the constraining power and degeneracies implied by 
the obtained constraints when we look at the final results, using the central values only as reference points. 

Below, we first detail the full parametrization of our likelihood analysis employed for the simulated YI and Y5 
analyses. 


With our prediction from Section O and mock data from Section HI 


B. Parameter Space 


The mock Y1/Y5 survey setup described in HIB yields a 20-dimensional parameter space. These parameters can 


largely be classified into cosmological, HOD, systematic effects, and growth scaling parameters. Here, we discuss how 
we set parameter defaults and priors for each parameter category, with references to relevant DES analyses on the 
SVAI data as well as external results serving as benchmarks. Table |l| lists the numerical values for parameter defaults 
and priors in detail. 

Cosmology - Eor cosmological parameters, we combine the Planck likelihood |31| with the likelihood that emerges 
from our pipeline, thereby enforcing Planck priors. Accordingly, our fiducial model takes Planck best fit parameters 
as defaults. With the addition of growth scaling parameters, our model apparently has three parameters (A^, and 
Ai) that shift the overall clustering strength for each lens bin. However, galaxy bias is not treated as a free parameter, 
but rather as a function of the halo and galaxy mass distribution (Equation III, and with the constraints on Ag from 


Planck, we are able to constrain Ai independently as initially suggested by Yoo and Seljak [5|. 

HOD - HOD priors represent the additional constraining power on HOD parameters that we expect to obtain 
from information not used by our current setup. Eor example, since our analysis does not use small-scale galaxy 
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Sector 

Parameter 

Fiducial Default 

Prior width (1-cr) 

Source 


Dm 

0.314 



Cosmology 

h 

0.673 

Planck likelihoods 

Planck Collaboration et al. [34] 


As 

2.15x10"® 




log Mmin 

12.36, 12.33 

0.09 


HOD 

log Ml 

<^log M 

13.69, 13.58 

0.32, 0.30 

0.05 

0.15 

Coupon et al. 1251 


a 

1.28, 1.37 

0.05 


Lens Photo-z 

^zL 

bzL 

0.02, 0.02 

0, 0 

O O 

d d 

Rozo et al. 1351 

Source Photo-z 

(^zS 

hzS 

0.08 

0 

0.01 

0.01 

Sanchez et al. [35], Bonnett et al. 1371 

Shear Calibration 

m-y 

0 

0.02 

Jarvis et al. [^, Clampitt et al. 1381 

Growth Scaling 

Ai 

1, 1 

Flat [0.5, 2.0] 

N/A 


TABLE I. List of parameters with their fiducial defaults and 1-a prior widths presented with the respective sources from which 
we draw these values. Entries with a pair of values represent parameters that vary between the two lens bins, while entries 
with a single value represent parameters that are global for both bins. Note that the prior widths are default settings for the 
conservative Y5 analysis; for analyses with different assumptions, subsets of parameter widths are varied as stated below. All 
mass values are units of Mq/H. 


clustering, we can imagine including HOD constraints from an independent small-scale galaxy clustering analysis as 
HOD priors. Or, if we later include small-scale galaxy clustering in our analysis, the expected strengthening of HOD 
constraints can be emulated by HOD priors in our current setup. Since an independent, HOD-focused analysis is yet 
to be performed on DES data, we use the results of the CFHTLS-Wide survey [5S] as a benchmark for the eventual 
DES HOD constraints. For fiducial defaults, we adopt the CFHTLS best-fit HOD parameters with a comparable 
luminosity and redshift selection. For priors, we consider two primary sets of assumptions. The first set, which we 
refer to as conservative, assumes that DES Y5 data will yield HOD constraints equivalent to the CFHTLS results, 
and use the CFHTLS 1-a uncertainties as default widths of Gaussian priors on the HOD parameters in the simulated 
Y5 analysis, as detailed in Table |Ij In the simulated Y1 analysis, we double the default prior widths for logMmin, 
log Ml, and a to reflect the relatively smaller sky coverage and shallower depth of the Y1 data stage. The second set, 
which we refer to as optimistic, assumes that DES Y1 and Y5 HOD constraints will scale with their increased sky 
coverages compared to CFHTLS, and uses HOD prior widths decreased by factors of (/sky,Yi//sky,CFHTLs) and 
(/sky,Y 5 //sky,CFHTLS ) for Y1 and Y5, respectively. These factors are roughly 2.7 and 6.1 for Y1 and Y5. One of 
the key issues in our analysis is how much information is needed about the HOD parameters in the quest to constrain 
the cosmological parameters and to study the effect of these priors on our eventual constraining power, we also 
carry out several conservative Y1 analyses where the the widths for Mmin, Mi and a are loosened by factors of 1.5, 
2.5, 3.5, and 5. 

Systematic Effects - By systematic effects parameters, we refer to the lens photo-z, source photo-z, and the 
multiplicative shear calibration parameters. We expect to understand the extent of photo-z errors present in DES 
catalogs from studies of spectroscopic subsamples and simulations, and this information can be incorporated into 
this analysis through photo-z priors. The lens photo-z modeling and priors adopted above are realistic for an LRG 
galaxy sample, and we anticipate the first application of our data to use DES redMaGiG [ 33 ] galaxies as the lens 
sample. The redMaGic galaxy sample is selected by fitting every galaxy to a red sequence template and establishing 
chi-squared cuts to enforce a constant comoving spatial density of galaxies over redshift, which by design allows for 
the selected galaxies to have tight and well-behaved (Gaussian) photo-z constraints. The “pessimistic” redMaGiC 
photo-z estimates are reported as a^i, = 0.015 and = 0 with ~ 1% catastrophic redshift failure rate, and we use 
conservative values of a^L = 0.02 and hzL = 0 as our defaults. For the photo-z precision of source galaxies, early 
photo-z results in DES data [ 35 ] suggests a^s = 0.08 and h^s = 0, which we use as defaults. We adopt Gaussian 
priors of width 0.01 for these four parameters, allowing both the bias and the variance of the photometric redshift 
estimates to be determined by the data, subject to modest priors on their ultimate values. In addition, we note that 
while tests of consistency between the lensing from different source redshift bins shows no discrepancies in the SVAl 
data (Clampitt et al. [38]), such tests do not account for an overall multiplicative bias that would affect all source 
bins equally. This multiplicative shear calibration parameter, rrij, is measured in Jarvis et al. |29j and found to be 
less than 2%. Thus, we assume = 0 as our fiducial default, and introduce a 0.02 (2%) Gaussian prior on this 
parameter. Finally, similar to the HOD priors, we carry out a mock Y1 analysis where the prior widths for systematic 





















12 


effects parameters are widened by a factor of 2.5 to gauge how they affect our final constraining power on Ai. This 
exercise is extended to a number of optimistic Y5 analyses with systematic effects prior widths of 0.5, 2, and 4 times 
the default width, as the optimistic Y5 scenario is expected to exhibit the strongest impact from systematic effects 
parameters. 

Galarxy abundance priors - Galaxy abundance, or the total number of galaxies in the survey, is calculated as 

Ng = ns j 

where fig is the solid angle subtended by the survey. The calculation of fig, from Equation has a different 
dependence on the mass function and the HOD than does the galaxy bias bg from Eq. so simply counting the 
number of galaxies in the survey will provide an additional constraint on HOD parameters. In particular, the number 
of galaxies in the survey breaks a problematic degeneracy between a and Ai. To implement galaxy abundance priors 
in our simulated examples, we will assume a generic 10% scatter in Ng for Y1 and 5% for Y5, and adopt corresponding 
Gaussian likelihoods into the analysis. Note that under the most ideal circumstances, there is only the Poisson and 
sample variance uncertainties on Ng. However, since galaxy selection is diluted by uncertainties in the photo-z and 
the mass-luminosity relation, we choose to adopt these conservative prior widths. 


C. Running and Verifying Chains 

To implement the MGMG, we use CosmoSIS [3S], a modular parameter estimation framework. Eor MCMC sampling, 
we make use of the emcee sampler m. an implementation of the affine-invariant MGMG ensemble sampler discussed 
in Goodman and Weare [41], using 190 walkers. Each ensemble iteration in our MGMG thus consists of 190 samples, 
one from each walker. In order to ensure that our chains are properly drawing independent samples from the likelihood 
space, we utilize measurements of the integrated autocorrelation time r as a criterion for testing convergence. Denoting 
the mean value of parameter pi in the t-th ensemble iteration as Pi(t), the autocorrelation function Ci(T) for that 
parameter with ensemble iteration lag T is given by 

C,{T) = HMt + T)- {p.))mt) - {p^))) ■ (22) 

The autocorrelation function is commonly normalized as 

p,(T)=C,(T)/C',(0), (23) 

which then yields the integrated autocorrelation time as 

(Tmax 

n = l+^MT). (24) 

T=1 

For a properly converged chain, reaches an asymptotic value and is stable with respect to Tmax, and this behavior 
then can be used as a heuristic signal for convergence. For our chains, stabilized values for different parameters 
range from 40 to 130 ensemble iterations. As suggested in [ID], we then consider the first few (around 10) ensemble 
iterations as burn-in, and choose to discard the first 1000 ensemble iterations. Parameter estimation is then performed 
on the following 900 ensemble iterations, consisting of 171,000 samples. In addition, we check that the acceptance 
fraction observed in our chain is stabilized to a reasonable value for the chosen region. 


V. RESULTS 

This section presents the results from the likelihood analyses described in Section [TV] revolving around “triangle” 
plots of 1-D and 2-D constraints of model parameters. In visualizing our results, we use a modified version of the 
triangle |42j Python package. In each subsection below, we will focus mostly on particular subsets of parameters. 

Large plots encompassing full parameter sets are pushed off to Appendix [^ In particular. Figs. [T^ and [Is] are our 
forecast parameter constraints from the default simulated Y1 and Y5 likelihood analyses. As discussed in Section 
|IV[ these results represent our conservative and optimistic estimates at the eventual DES constraining power on the 
growth function for the respective data stages. The takeaway is that parameter constraints are very well centered 
with respect to their true values, an indication that the 20-parameter MGMG is working well. Note that we are not 
showing constraints on the cosmological parameters, as these are largely dominated by Planck priors. Therefore, we 
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suppress these columns but come away with the knowledge that for the cosmological parameters that we consider - 
O.M, h, and Ag - we expect CMB constraints to be dominant over constraints from combining small scale leasing and 
large scale clustering. Also, as the two lens bins show very similar parameter behaviors, we only show contours for 
the first lens bin and simply tabulate results for the second lens bin. 

In all triangle plots shown below, the panels along the diagonal correspond to the 1-D probability distributions for 
each parameter with dotted vertical lines at the 16th and 84th percentiles, while the off-diagonal panels show the 2-D 
Gaussian 1-a confidence contours for the corresponding pair of parameters. The light blue lines and squares represent 
the fiducial default parameter values used to generate the simulated measurement vector, i.e. the “true” parameter 
values. Numerical values for marginalized l-cr bounds are listed in Table [III 


A. Default Conservative and Optimistic Results 

Let us begin with results under the two default - conservative and optimistic - assumptions for the simulated Y1 
analysis. In Fig. we present our forecast parameter constraints on the HOD and the growth scaling parameters 
from those two analyses. A key issue for our study is the extent of the correlation between Ai, representing the 
amplitude of matter fluctuations in the lens bins, and the HOD parameters. If there were no degeneracy, then the 
analysis could be carried out without any dependence on HOD modeling. Fig. shows that this is not the case, i.e. 
that the HOD parameters are correlated with A^. In particular, Aj are quite degenerate with the two parameters that 
quantify the satellite galaxy abundance. Mi and a. 

The difference between the two analyses is in widths of HOD priors, where the conservative analysis assumes 
widths twice as large as CFHTLS constraints, while the optimistic analysis assumes widths smaller by a factor of 
(/sky,Yi//sky,CFHTLs)^^^ (roughly 2.7), resulting in a net difference in HOD widths by a factor of 5.4. For central 
HOD parameters Mmin and CTiogM: we observe that the difference between the conservative and optimistic constraints 
is smaller than the difference in the prior widths. For satellite HOD parameters Mi and a, we observe the difference in 
constraints closely following the difference in the prior widths, indicating that these constraints are largely prior-driven. 
From the two analyses, we project 7.9% (conservative) and 4.7% (optimistic) l-cr error bars on Ai. 



FIG. 4. Forecast constraints on the HOD and growth scaling parameters from the conservative (blue) and optimistic (red) 
Y1 analyses. Diagonal blocks represent marginalized 1-D parameter constraints, and off-diagonal blocks represent 2-D 1-a 
confidence ellipses for corresponding parameters. 


In Fig. we compare the results from the conservative and optimistic Y5 analyses. The conservative analysis 
assumes widths equivalent to CFHTLS constraints, while the optimistic analysis assumes widths smaller by a factor 
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1 

of (/sky.Ys/Zsky.CFHTLs) ^ Or roughly 6.1. We observe similar parameter behaviors as in the Y1 counterparts, and 
project 3.9% (conservative) and 2.3% (optimistic) l-cr error bars on Ai. 



FIG. 5. Forecast constraints on the HOD and growth scaling parameters from the conservative (bine) and optimistic (red) 
Y5 analyses. Diagonal blocks represent marginalized 1-D parameter constraints, and off-diagonal blocks represent 2-D 1-a 
confidence ellipses for corresponding parameters. 

Let us now turn our attention to the systematic effects parameters. In Fig. we present our forecast parameter 
constraints on the systematic effects and the growth scaling parameters from the conservative and optimistic Y1 
analyses. As opposed to the HOD parameters. Fig. [^shows that there are no notable degeneracies between systematic 
effects parameters and our parameter of interest Ai. In Fig. we show the conservative and optimistic Y5 constraints 
for systematic effects and growth scaling parameters, and observe identical parameter degeneracies. Also, note that 
for both results how small the difference in constraints is between conservative and optimistic results. This implies 
that changing HOD priors has only a marginal effect on systematic effects constraints, i.e. that HOD and systematic 
effects parameters show little degeneracy between them. This is also observable in Figs. [T^and[T5| 

To further explore the effect of degeneracies between different parameter subsets and Ai, as well as to gauge how 
much constraining power is coming from the analysis itself as opposed to assumed priors, we perform a series of 
exercises in relaxing the priors for different parameter subsets. 


B. Relaxing HOD Priors 

To get a sense of how important the priors on the HOD parameters are, we perform a “relaxed HOD priors” 
simulated Y1 analysis, where the width of HOD priors are widened by a factor of 2.5 compared to the conservative 
default widths. In Figure]^ we compare the conservative Y1 constraints for HOD and growth scaling parameters, 
as presented in Figure]^ against Y1 constraints with relaxed HOD priors. Our constraining power on Mmin is only 
mildly degraded despite the relaxed priors, implying that our analysis constrains largely by itself without relying 

on external priors. However, Mi and a are not as well constrained by the new, wider priors, implying that the priors 
are driving our constraining power on these parameters. Since these are the parameters that are most degenerate 
with Ai, it is not surprising that the constraint on Ai loosens by a factor of 2, with the 7.9% error bar on Ai using 
nominal priors degraded to 16.1% with the looser priors. 

In Fig. we go further to quantify the extent to which Mi and a impact Ai. we compare our final constraints 
on the growth scaling parameters Ai ^2 from a number of different Y1 analyses as a function of the prior width on 
the HOD parameters. In addition to the default conservative and 2.5x widths, we consider prior widths equal to 
0.5, 1.5, 3.5, and 5 times the default conservative width. The minimum mass M„iin is relatively well constrained 
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FIG. 6. Forecast constraints on the systematic effects and growth scaling parameters from the conservative (blue) and optimistic 
(red) Y1 analyses. Diagonal blocks represent marginalized 1-D parameter constraints, and off-diagonal blocks represent 2-D 
1-a confidence ellipses for corresponding parameters. 


by the data regardless of priors, so the main effect of varying HOD priors is on the satellite HOD parameters. The 
result we observe, as presented in Fig. is a linear relationship between the prior widths and the constraints on 
Hi 2 , which confirms the strong degeneracy between satellite HOD parameters and Hi 2 . The lessons from Fig. [^are 
straightforward: the constraints on the cosmological parameters of interest will be limited by our ability to constrain 
the satellite HOD parameters using other measurements. 


C. Relaxing Systematic Effects Priors 

It is also important to understand the effect of systematic effects - photometric redshift uncertainties and biases for 
both source and lens galaxies and multiplicative shear calibration - on the cosmological constraints. This understand¬ 
ing could be used to implement scientific requirements on shear for joint analyses of these types (they may be looser 
than those needed for cosmic shear) and to estimate the number of spectroscopic redshifts needed to reduce photo¬ 
metric redshift errors. To gauge the impact of systematic effects priors on our final constraining power, we carry out 
another conservative Y1 analysis with relaxed systematic effects priors, similar to the relaxed HOD priors case above, 
where we widen the widths of systematic effects priors by a factor of 2.5 compared to the default conservative widths. 
In Figure [TOl we compare the default conservative Y1 results, as presented in Figure]^ against the Y1 results with 
relaxed systematic effects priors. Perhaps the main takeaway is the bottom right panel, which shows that constraints 
on Hi are degraded minimally (from 7.9% to 9.4%) in this case, when the prior constraints on systematic effects 
parameters are significantly relaxed. This modest degradation is due partly to lack of degeneracy between Hi and 
most of the systematic effects parameters; note for example the flattened ellipses in the bottom row in columns for cr^g 
and bzL. Even though these nuisance parameters are not well-constrained by the data (blue ellipses are much wider 
than red), they are not degenerate with Hi, so they have a limited effect on the final growth constraints. By contrast, 
the data does constrain (Jzl, the scatter in the lens photometric estimates, quite well even without an external prior. 

While the systematic effects parameters had only a small effect in the eventual growth constraints for the conserva¬ 
tive Y1 analysis, we anticipate that its relative impact will be bigger in the Y5 scenario, especially for the optimistic 
Y5 analysis where we employ the tightest HOD priors. To test the significance of systematic effects in the Y5 case, 
we perform an exercise similar to that for the HOD priors and compare our optimistic Y5 constraints on the growth 
function with varying prior widths for systematic effects parameters. In Fig. El we present the fractional uncertainty 
on the growth scaling parameters Hi _2 with systematic effects priors 0.5, 1, 2, and 4 times the default width. We 
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FIG. 7. Forecast constraints on the systematic effects and growth scaling parameters from the conservative (blue) and optimistic 
(red) Y5 analyses. Diagonal blocks represent marginalized 1-D parameter constraints, and off-diagonal blocks represent 2-D 
1-a confidence ellipses for corresponding parameters. 


observe that even under the optimistic Y5 scenario the degradation from relaxing systematic effects priors to growth 
constraints is weak. As we relax the systematic effects priors, we also found that the lens photo-z parameters (ct^l 
and bzh) and the shear calibration parameter m-y constraints do not weaken as much as the relaxation in priors, while 
the source photo-z parameters {a^s and bzs) exhibit changes in constraints that closely follow the relaxed priors. In 
terms of impact of priors, this indicates that the pipeline is constraining the lens photo-z and the shear calibration 
parameters by itself, while it is relying on priors to constrain the source photo-z parameters. In terms of growth 
constraints, this indicates that those same parameters that we are constraining without heavily relying on priors, 
i.e. lens photo-z and shear calibration parameters, are dominant over those that we constrain largely by priors, i.e. 
source photo-z parameters, in their impact on our final constraints on the growth function. Finally, the decreasing 
trend that we observe in Fig. 11 as we tighten our systematic effects priors beyond the default width implies that 


better understanding and constraining of systematic effects will be important in achieving the tightest possible growth 
constraints from DES. 
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FIG. 8. Forecast constraints on the HOD and growth scaling parameters from the default conservative (red) and relaxed HOD 
priors (blue) Y1 analyses. The relaxed HOD priors are a factor of 2.5 weaker than the default priors. 
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FIG. 9. Fractional l-cr uncertainty on Ai for different Y1 analyses, plotted against HOD prior widths as a factor of the default 
conservative widths. 
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FIG. 10. Forecast constraints on the systematic effects and growth scaling parameters from the default conservative (red) and 
relaxed systematic effects priors (blue) Y1 analyses. The relaxed systematic effects priors are a factor of 2.5 weaker than the 
default priors. 
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FIG. 11. Fractional 1-cr uncertainty on Ai for different simulated Y5 analyses, plotted against systematic effects prior widths 
as a factor of the default widths. 
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VI. DISCUSSION 

In this paper we demonstrate an implementation of the joint-analysis pipeline for combining galaxy clustering and 
galaxy-galaxy leasing measurements from photometric surveys. In preparation for DES data analyses, our modeling 
includes the expected key systematic effects of photometric redshift estimates, shear calibration, and the galaxy- 
luminosity mass relationship, covering a 20-dimensional parameter space. We show that a joint analysis of large-scale 
w{9) and small-scale 7t(0) can conservatively/optimistically constrain the growth function D{z) to within 7.9%/4.8% 
with DES Y1 data and to within 3.9%/2.3% with DES Y5 data across two different redshift bins of 0.3 < z < 0.4 
and 0.4 < z < 0.5. These forecasts can be put in the context of existing constraints on the growth function using 
the abundance of galaxy clusters, weak leasing shear correlations, and redshift space distortions in galaxy clustering. 
Some recent results include: 

• Galaxy Clusters [ 5 ]: erg = 0.83 ± 0.04 

• Weak Tensing gS]: 0-8(11^/0.27)“ = 0.774/°;°^^, a = 0.46 ± 0.02 

• Weak Tensing [44]: o-8(IIm/0.3)° ® = 0.81 ± 0.06 

• Redshift Space Distortions gS]: o-s = 0.76 ± O.ll 

• Redshift Space Distortions gSj: /(z)o-8(z) = 0.413 ± 0.080,0.390 ± 0.063,0.437 ± 0.072 at z = 0.44,0.6,0.73 

• Redshift Space Distortions g7]: fizes)crs{zes) = 0.419 ± 0.044 at z^s = 0.57 

Overall, current constraints are at the 10% level, implying that DES will soon produce cutting-edge constraints on 
the growth of structure in the universe with tomography reaching far back in cosmic time, even under conservative 
assumptions. Fig. [T^ shows the bounds on the ACDM growth function obtained for the two lens bins under the DES 
Y1 and Y5 specifications considered in this analysis. 



FIG. 12. Forecasts for the DES 1-cr bounds on the growth function D{z) in two different redshift bins, presented for different 
assumptions on data stage and parameter priors. Points are offset for visibility. 

An important conclusion is that the HOD parameters are degenerate with our parameters of interest, i.e. the 
growth scaling parameters A,, but the systematic effects parameters are at most weakly degenerate with A^. By 
comparing results drawn under different prior settings, we conclude that the final constraining power on the growth 
function will be driven by our ability to constrain HOD parameters, especially the satellite HOD parameters, as these 
parameters are both strongly degenerate with A^ and relatively unconstrained without priors. On the other hand, 
we observe that the central HOD parameters are well constrained without contribution from priors, and also that 
the systematic effects parameters are either well-constrained or only weakly affecting the final constraining power on 
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Ai- The default results are thus strongly driven by the priors on the satellite HOD parameters, but as these prior 
constraints encompass both conservative and optimistic estimates of the eventual constraining power coming from 
DES, we believe looking at both the conservative and optimistic forecasts yields reasonable estimates for the ultimate 
results we can expect from analyses of DES Y1 and Y5 datasets. 

The HOD degeneracies also illustrate limitations in the two-step analysis proposed in Yoo and Seljak m of (1) 
determining the (largely cosmology-independent) mean halo mass of the galaxy sample from stacked small-scale 
leasing measurements and (2) analyzing large-scale galaxy clustering using galaxy bias inferred from the obtained 
mean halo mass to determine the amplitude of the underlying matter clustering. In this method, the connection from 
the first to the second step, and consequently the determination of galaxy bias, hinges on a single representative value 
- the mean halo mass obtained from galaxy-galaxy leasing. Thus, if galaxy samples with similar mean halo masses can 
exhibit varying galaxy biases, the two-step approach becomes sub-optimal. And the HOD degeneracies suggest that 
such a situation is entirely possible. In Fig. we show 10,000 random HOD configurations with mean halo masses 
within 1% of our fiducial default HOD for lens bin 1, along with the derived galaxy bias from those configurations. 
Note that these are not results from an MCMC analysis, but simply random HOD’s within a relatively narrow range 
of parameter values that yield the desired mean halo masses. Even with this extremely tight requirement in mean 
halo mass, different galaxy samples exhibit a much wider scatter (up to 10%) in galaxy bias, as shown in the panel in 
bottom right corner. This result suggests that the approach employed in our analysis, i.e. a consistent HOD modeling 
of a given galaxy sample that propagates to predictions for both galaxy clustering and galaxy-galaxy leasing, is a 
more optimal form of combined probed analysis. 



FIG. 13. Plot of 10,000 random HOD configurations with mean halo masses within 1% of the fiducial default value for lens bin 
1, presented with the galaxy biases derived from these configurations. The plot ranges are set to match the range of parameter 
values used to generate random HOD configurations. Light blue squares and lines mark the fiducial default HOD for lens bin 
1, as presented in Table |I] 


Based on the lessons learned from this study, our current implementation will undergo a number of key improvements 
in the near future. The most salient improvement will be incorporating small-scale galaxy clustering information. 
Small-scale galaxy clustering is highly sensitive to satellite galaxies, and thus will allow for tight constraints on the 
satellite HOD parameters. With this improvement, all of our HOD constraints will be data-driven, and our analysis 
will be self-sufficient without relying on HOD priors. In addition, we expect further validation of our assumed HOD 
model from a separate DES analysis on HOD modeling. This improved pipeline is anticipated to analyze the DES 
SVAl and Y1 data and produce interesting constraints in the near future. 
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Appendix A: Forecast Figures 



FIG. 14. Marginalized ID parameter constraints with dotted vertical lines at ±la (diagonal) and 2D 1-a confidence ellipses 
(off-diagonal), representing the conservative (blue) and optimistic (red) Y1 parameter constraint forecasts for the first lens bin. 
Vertical and horizontal axes represent the 4 HOD parameters, 5 systematic effects parameters, and the growth scaling parameter, 
respectively. Light blue lines and squares correspond to the true values used in generating the simulated measurement vector. 
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FIG. 15. Marginalized ID parameter constraints with dotted vertical lines at ±1 ct (diagonal) and 2D 1-cr confidence ellipses 
(off-diagonal), representing the conservative (blue) and optimistic (red) Y5 parameter constraint forecasts for the first lens bin. 
Vertical and horizontal axes represent the 4 HOD parameters, 5 systematic effects parameters, and the growth scaling parameter, 
respectively. Light blue lines and squares correspond to the true values used in generating the simulated measurement vector. 























































































